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Abstract 

We outline a novel linear instability that may arise in the dead-zones of protostellar disks, and possibly 
the fluid interiors of planets and protoplanets. In essence it is an axisymmetric buoyancy instability, but one 
that would not be present in a purely hydrodynamical gas. The necessary ingredients for growth include a 
negative radial entropy gradient (of any magnitude), weak magnetic fields, and efficient resistive diffusion 
(in comparison with thermal diffusion). The character of the instability is local, axisymmetric, and double- 
diffusive, and it attacks lengths much shorter than the resistive scale. Like the axisymmetric convective 
instability, it draws its energy from the negative radial entropy gradient; but by utilising the diffusing 
magnetic field, it can negate the stabilising influence of rotation. Its nonlinear saturated state, while not 
transporting appreciable angular momentum, could drive radial and vertical mixing, which may influence 
the temperature structure of the disk, dust dynamics and, potentially, planet formation. 
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1. Introduction 

In this paper, we consider the dynamical behavior of a gaseous system characterized by rotation (with a 
strong positive angular momentum gradient), a weak magnetic field, and a small negative entropy gradient. 
As is well known, these ingredients lead to the magnetorotational instability (MRI), possibly in its more 
general convective form (Balbus & Hawley 1991, Balbus 1995). It might appear that when resistivity 
is present only the longest wavelengths can remain unstable. If the resistivity is too large then these 
wavelengths can be longer than the size of the system, and global stability results. (Note that the small 
negative entropy gradient is stabilized in the hydrodynamical limit by strong rotation.) 

Surprisingly, this is not an accurate description of the dynamical behavior of this system. Here we show 
that the presence of resistivity activates unstable buoyancy modes that rely on the arbitrarily small adverse 
entropy gradient as their free energy source. Moreover, these unstable modes lie on the same branch of the 
dispersion relation as the MRI, to which they revert at small wavenumbers. What is particularly striking is 
that at very large wavenumbers, resistivity never damps these modes. Ultimately, when thermal diffusion is 
considered, there is in fact a critical wavenumber above which damping occurs. But, at least in protostellar 
disks and planetary interiors, this still allows for a very broad range of unstable large wavenumber modes. 
(The presence of both Ohmic resistivity and a much smaller thermal diffusivity is what gives this instability 
its 'double diffusive' character.) 

In this introductory paper, our interest is in a formal presentation of the properties of the instability in 
the simplest possible setting. Our fiducial model will be a protostellar disk. Such systems are both cold and 
dense; excepting the hot regions near the central star and in the surface layers, much of the gas is poorly 
ionised, and, as a consequence, the electrical conductivity too low to develop classical MRI turbulence (Blaes 
and Balbus 1994, Gammie 1996, Igea and Glassgold 1999, Sano et al. 2000, Fleming and Stone 2003, Ilgner 
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and Nelson 2006, Wardle 2007, Okuzumi 2009, Turner and Drake 2009). Most protostellar disk models 
hence consist of an 'active' well-ionized envelope, in which the MRI could in principle be present, encasing 
an extensive body of poorly- ionised gas which remains 'dead', at least as far as the MRI is concerned. The 
dead region nevertheless may support other dynamics, such as spiral density waves excited by the turbulence 
in the surface layers (Fleming and Stone 2003), streaming instability caused by settling dust grains (Youdin 
and Goodman 2005), and flows driven by torques applied by large-scale magnetic fields (Turner and Sano 
2008). Generally, however, dead-zones have been considered magnetically 'deactivated'. Nevertheless, even 
in the very resistive heart of the dead-zone we will see that magnetic fields can trigger local small scale 
instability. 

The analysis will be explored in a local model, the details of which we present in Section 2. In Section 
3, we investigate the dynamics of linear axisymmetric Boussinesq disturbances by solving their fifth-order 
dispersion relation numerically and by obtaining asymptotic estimates of the growth rates in the double- 
diffusive limit. The instability is subsequently related to the idea of magnetostrophic balance. In Section 4 
our conclusions are drawn and future work outlined. 



2. Physical and mathematical outline 

Consider a weakly ionized region of a protostellar disk, a nominal 'dead-zone'. A negative radial entropy 
gradient is assumed to be present, but the details of how this gradient is established and maintained do not 
concern us directly in this calculation: they may plausibly be present due to a combination of outwardly 
decreasing heating from the central star and an outwardly increasing opacity. The entropy gradient is 
presumed to be small relative to the angular momentum gradient of the disk. More specifically, the radial 
Brunt- Vaisalla (BV) frequency 

2 1 dP d In 
Nr ~ 1P dR OR W 
is assumed to be small compared with the epicyclic frequency 

Here, R is the radial coordinate in a standard cylindrical system (R,(f),z), p is the mass density, P is the 
gas pressure, VL(R) is the angular velocity, and 7 is the adiabatic index (ratio of constant pressure specific 
heat to constant volume specific heat). Recall that k — in a Keplerian disk, so the magnitude of Nr is 
assumed much less than an orbital frequency (Nr itself is, of course, imaginary here). 

We consider the stability of axisymmetric disturbances. If we were to treat this highly resistive system 
as purely hydrodynamic, we would conclude that it is stable: there is no instability if N R + n 2 > (the 
Hoiland-Solbcrg criterion), which for small \Nr\ < Q is certainly true. On the other hand, MHD processes 
are generally also stabilized when the resistivity is high. The classical MRI, for example, is banished to 
unfeasibly long wavelengths in the presence of strong Ohmic dissipation. Consequently, as noted in the 
Introduction, one might be tempted to conclude that adding a very small negative entropy gradient to a 
magnetised resistive disk would make very little difference to the axisymmetric stability. In fact, it makes a 
dramatic difference. 

Consider two fluid elements at different vertical locations embedded in a negative radial entropy gradient. 
A weak magnetic field tethers these two elements to one another. When these elements are radially displaced, 
they will try to bring the magnetic field with them, but the field diffuses because of resistivity (which would 
quell the classical MRI). Nevertheless, because of the transient magnetic torque, some finite amount of 
angular momentum will have been exchanged between the two elements in the process. As a consequence, 
each element will settle into a new orbit at a new radial location — and, most critically, in a new entropy 
environment. If thermal diffusion is too slow to equilibrate the temperature of the displaced fluid with 
its surroundings, the elements will continue to move radially, because they are now buoyantly unstable. 
Subsequently, the diffusing field will ensure that these elements will magnetically connect to new fluid 
parcels, and further angular momentum will be exchanged. Thus the cycle continues and an instability is 
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at hand: an instability that depends upon the transfer of angular momentum between magnetized fluid 
elements, yet whose seat of free energy is not the shear, but rather the unstable thermal stratification. Note 
that this instability is double-diffusive, relying on the speedy diffusion of angular momentum (accomplished 
by diffusing magnetic field) and the slow (or negligible) diffusion of heat. This field diffusion counters 
the stabilising tendency of an angular momentum gradient by breaking the constraint of specific angular 
momentum conservation. 

Though buoyantly driven, the instability clearly has features in common with the MRI. On sufficiently 
small vertical wavenumbers (sufficiently spaced fluid elements) the influence of resistive diffusion will be 
minor and we recover the MRI. In fact, we will show that the double-diffusive instability and the MRI are, 
in effect, two sides of the same instability. In the presence of a negative radial entropy gradient, one of the 
modes that emerges is, at small wavenumbers, the classical MRI and, at large wavenumbers, the double 
diffusive instability. The instability is suppressed only at extremely large wavenumbers, when thermal 
conduction becomes important. Thus the MRI and the resistive double-diffusive instability are intimately 
linked. 

In the following sections this idea is explored using a number of approaches, mainly within the framework 
of the linearised Boussinesq equations of resistive MHD. (An important omission, which we will examine in 
a later study, is the electromotive force associated with the Hall term in the induction equation.) The local 
modes under consideration satisfy a fifth-order dispersion relation, which we will study in some detail. It is 
possible to obtain a number of clean numerical and analytical results from which further physical insights 
emerge. 

2.1. Model equations 

Our model employs the resistive MHD equations, which comprise the continuity, momentum, and entropy 
equations, and the magnetic induction equation with Ohmic diffusion. As noted, the Hall effect is also of 
practical importance, but it is omitted in this first analysis. Ohmic resistivity will exceed viscosity and 
thermal diffusion by several orders of magnitude, hence the viscous stress is dropped. At least initially, 
however, it is helpful to retain the thermal diffusion term, which we model as the divergence of a radiative 
energy flux. The equations are 

^ + v • Vp = -pV • v, (3) 
<9v „ l„/ n B 2 \ ^ (B-V)B 

9B 

— + vVB = B-Vv-BV-v + ?) V 2 B, (5) 

E +V ' V(T ) = ~ V ' F + ^ V x B l 2 ' ( 6 ) 

Here v is the velocity, E the internal energy density, $ the gravitational potential, B the magnetic field, 
and a = InPp -7 is the entropy function. The resistivity is represented by 77, and the radiative energy flux 
density by F which is defined through 

where T is temperature, n op is opacity, and 03 is the Stefan-Boltzmann constant (the latter two are not to 
be confused with the epicyclic frequency k and the entropy function a). Finally, the equation of state is of 
an ideal gas which gives 

E=^—P. (8) 

In this study, we restrict the analysis to current-free initial equilibria, i.e. V x B = 0, and the last term of 
equation (6) will not be used. 
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2.2. Linearised equations 

It is assumed that the disk is in near Keplerian equilibrium, with a small pressure gradient supplementing 
the star's gravity in the radial force balance. The equilibrium is v = RCl(R)e^,, p = p(R), P = P(R) 
and S = S(R). It is also assumed from the outset that the radial gradients of P and S are negative and 
relatively small. In addition, weak magnetic fields thread the disk, but they are too small to influence the 
equilibrium balances. For simplicity, the magnetic field takes a uniform vertical configuration: B = _B e z . 

In order to examine local disturbances, we adopt the shearing sheet model (Goldreich and Lynden-Bcll 
1965), which represents a small block of disk as a Cartesian box with (x, y, z) denoting the radial, azimuthal, 
and vertical coordinates. Inside the box p, dP/dR, and dS/dR are assumed to be constant. 

This equilibrium is perturbed by a planar axisymmetric Boussinesq disturbance, p', v' x , v' y , B' x , B' y , 
proportional to e tkz+st , the linearised equations of which read 

i „^ i k 2 , p' ( dP\ ikBo^, 

. v< = 2n„ „<„ e, - ^ K e» + I ) ^ e, + — B', (9) 



sB ' = ifcB ^ + (^)„ i!;e '-'" t2B ' 



(10) 



where the subscript indicates evaluation at the point on which the shearing sheet is anchored. Hereafter, 
the subscript will be dropped. The thermal diffusivity £ has been introduced which is defined by 

^ = 3 7 K op pP (12) 
Notice that in the entropy equation (11), we have used the relation 

r = -(T/p)p\ (is) 

which comes from the perturbed ideal gas law. The pressure perturbation is identically zero for our pla- 
nar disturbances. But, in any case, Boussinesq perturbations are in near pressure equilibrium with their 
surroundings, and so the fractional pressure perturbation is negligible in the ideal gas law. 

3. Analysis 

Equations (9)-(ll) present an algebraic eigenvalue problem for s, which yields the following dispersion 
relation 

s 5 + a 4 s 4 + a 3 s 3 + a 2 s 2 + a\ s + a = 0, (14) 

with the coefficients 

a 4 = 2r,k 2 + £/c 2 , (15) 
a 3 = N 2 l + K 2 + 2v 2 A k 2 +r 1 2 k i + 2r]£ fc 4 , (16) 
a 2 = 2N 2 R nk 2 + 2(r, + Z)v 2 A k 4 + (2r? + £)k 2 k 2 + r? 2 £ fc 6 , (17) 

ai = v\ k 2 +v\k 4 + (n 2 k 2 + v\) k 2 N 2 R + (2^ + rj 2 ) fc 4 k 2 + 2^v\ fc 6 , (18) 

a = N 2 R r, v\ fc 4 + i v\ fc 6 + ft fc 6 n 2 + tf A fc 4 (J^j , (19) 

where the Alfven speed is defined through va = B n /y/4np. Equation (14) is a special case of the dispersion 
relation calculated by Menou et al. (2004), which also accounts for viscosity, vertical variation in ft, P, and 
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S, and radial wavenumbers. On the other hand, by taking r\ = £ = we recover the ideal MHD relation of 
Balbus (1995), which yields the two convective modes and the two convective-MRI modes. The dispersion 
relation in this case is fourth order, which means that the inclusion of magnetic diffusion (and thermal 
diffusion) leads to the emergence of a new fifth mode. 

From simply inspecting the sign of the last coefficient ao it is already clear that the system possesses 
novel stability behaviour when £ is small. If we consider scales upon which the thermal diffusion has little 
influence, the stability discriminant is the first term on the right side of (19). It follows that stability is 
linked to the sign of N^, which means that the Schwarzchild criterion is recovered — the stabilising effect 
of rotation has vanished. Notice also that the discriminant combines the three elements essential to the new 
instability: resistivity, entropy stratification, and magnetic fields. If any one of these is missing a vanishes 
and the systems' stability properties revert to the familiar situation of Balbus (1995), which is governed by 
the coefficient a±. 

In the following subsections we adopt a number of approaches in order to describe the mathematics and 
physics of this double-diffusive mode more fully. First, we briefly discuss the main parameters that appear 
in the model and estimate realistic values they may take in protostellar disks. Next we solve the dispersion 
relation Eq. (14) numerically, and obtain an analytical asymptotic result in the regime of interest. Following 
that, we show how the mode relies on a generalisation of 'magnctostrophic balance' in order to function. 
Finally, we give a more physical explanation of the mechanism of instability. 

3.1. Fiducial parameters and lengthscales 

Time and wavenumber units are chosen so that ft = 1 and the Alfven wavenumber, defined through 
kA = FI/va, is also 1. The Alfven length is denoted by I a and is defined by I a = 2n/kA- Consequently, 
the governing parameters of the dispersion relation comprise the scaled squared BV frequency iV^/fi 2 , the 
scaled squared epicyclic frequency k 2 /Q 2 , the Roberts number q, and the Elsasscr number A. The latter 
two are defined by 

Q=~, A =4- (20) 
n nil 

In order to estimate q we use the following formulas for the resistivity and thermal diffusivity 

r? = 2.34xl0 17 f-I—) 1/2 (^—] cm 2 s-\ (21) 



100 KJ V X e 



where x e denotes electron fraction (Blaes and Balbus 1994). In the latter expression we have used 7 = 7/5, 
and endowed the gas with a mean molecular weight of 2.3. According to the commonly used minimum mass 
model, at a few AU the disk may be characterised by T ~ 100K and p <~ 10~ 9 gcm~ 3 (Hayashi et al. 1985), 
while at these low temperatures K op ~ lcm 2 g _1 (Henning and Stognienko 1996). The ionisation fraction 
x e , on the other hand, is poorly constrained and extremely sensitive to dust grain size and abundance, the 
sources of ionisation, and vertical location. Various model show it ranging over many orders of magnitude: 
from lO" 13 to values as low as 10~ 19 (Sano et al. 2000, Ilgner and Nelson 2006, Wardle 2007, Turner and 
Drake 2009). Given this uncertainty, we estimate the Roberts number by q ~ a: e /10~ 9 . This suggests an 
upper limit for q of 10 -4 and an extreme lower limit of 10~ 10 . 

The Elsasser number requires knowledge of the strength of the latent magnetic field, which we parame- 
terise by the plasma beta: 

2 c 2 

13 = , (23) 

where c is the gas sound speed. This gives 



(24) 



with H — c/Vl the disk scale height. The bracketed term in Eq. (24) is the magnetic Reynolds number, for 
which we can obtain a rough estimate using equation Eq. (21) and assuming H = OAR and R w 5 AU. For 
T ~ 100K, this presents A <~ {3~ 1 (x e /10~ 15 ). So one may expect A to vary between 10" 1 and 10~ 6 or even 
lower. 

Finally, it is necessary to estimate the squared BV frequency iVj|, which represents the entropy gradient 
from which the double-diffusive modes derive their energy. Most a-disk models in fact yield a positive BV 
frequency (Balbus and Hawley 1998), but this is because they omit the central star's radiation field which 
will tend to drive a negative radial temperature gradient. However, it is difficult to quantify how this central 
heating is mediated by the turbulent inner radii of the disk, and the dead-zone itself, and so it is unclear 
which estimates one should use for the magnitude of N^. We simply assume that Aj^/fi 2 takes values 
between —0.1 and —0.01. Because of the pressure gradient the background flow fi is not strictly Keplerian, 
however the deviation should be minimal and so we set n w fi. 

The Alfven length scales as 

Therefore modes that possess a wavelength comparable or smaller than the Alfven length may fit comfortably 
into the disk when the fields are sufficiently sub-equipartition. In contrast the resistive length, defined by 
l v = i]/va, scales like 

l v <~ A^ 1 /^. 

There can be no classical-MRI on scales less than l n . Using the above estimate for A we derive l v <~ 
/3 1 / 2 (10~ 15 /x e ) H. Then, for realistic parameter values (3 > 10 2 and x e < 10 -14 , the following ordering is 
established: 

A (MRI) > l n > H > l A . (25) 

Thus the classical MRI modes will not occur in the main body of the disk in most cases, giving rise to a 
dead-zone. The double-diffusive instability, on the other hand, will be present on scales from l n down to a 
critical length l c upon which radiative diffusion becomes important. This critical length can be estimated 
by setting the constant coefficient ao in (19) to zero. We find 

Unless magnetic fields and the entropy gradient are exceptionally weak, this scale l c should be less than the 
scale height. We conclude that unstable double-diffusive modes will generally fit into the dead-zone. 

3.2. Numerical solutions 

In this section, the dispersion relation (14) is solved numerically. Once the parameters A, K 2 /f2 2 , N^/fl 2 , 
and q are stipulated, we may compute the growth rate s as a function of vertical wavenumber k. Throughout 
we set q to be a constant, q = 10~ 6 , and Ac/fi = 1, while varying A and A|./f2 2 . Hence, the Elsasser number 
A may be interpreted as a surrogate for magnetic field strength. The results of this subsection are limited to 
the MRI/doublc-diffusivc mode, which is the only mode that grows, and we neglect the other four decaying 
modes. 

Figure 1 illustrates the remarkable stability behaviour unlocked by the combination of resistivity and a 
negative entropy gradient. The purple dotted-dashed curve represents the MRI in ideal MHD (r] = but 
with Nft ^ 0) and the red dashed curve represents the MRI in resistive MHD but with no radial stratification 
(77 7^ and Nft = 0). In both cases instability is extinguished on relatively long scales: near an Alfven length 
in the former case, and near the resistive length in the latter case. (Because here A = 1 the two scales are 
comparable.) However, when resistivity and stratification are present (77 7^ and 7^ 0) instability 
extends to much shorter scales, as witnessed by the blue solid curve in Fig. 1. We identify this 'extension' 
of the MRI as the double-diffusive instability. The MRI transforms smoothly into the new instability upon 
a transitional scale comparable to l n ; both instabilities occur on the same branch of the dispersion relation. 
When the MRI mechanism is quenched (because of resistivity), the double-diffusive mechanism takes its 



6 




Figure 1: Growth rate s of the unstable MRI/double-diffusivc mode as a function of k for three illustrative cases. The purple 
dashed-dotted line represents the ideal MHD configuration A = oo with N^/Q. 2 = —0.1. Here the unstable mode is the 
convective-MRI and it is extinguished when k/kj± > \/3- The red dashed line represents resistive MHD in the absence of a 
radial entropy gradient; here A = 1 and = 0. Instability is extinguished on wavelengths roughly below the resistive scale. 
The solid blue line represents the case with both resistivity and the negative entropy gradient: A = 1 and N^/D. 2 = —0.1. As 
is plain, instability continues into short scales; this is the double- diffusive instability. In the non-ideal case, the ratio of thermal 
to resistive diffusivity is set to q = 10~ 6 . 
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Figure 2: Growth rates of the double-diffusive mode s as a function of k for different Elsasser number A. The blue line represents 
the A = 0.1 case, the purple line, A = 0.05, and the red line, A = 0.01. The squared BV frequency is TV^/O 2 = —0.1 for these 
cases and the Roberts number is q = 10~ 6 . At very long scales, k -C k v ~ Afc^, we recover the familiar MRI mode, but as k 
approaches k v this mode morphs into the double-diffusive mode which persists until k c . 
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Figure 3: Growth rates of the double-diffusive mode s as a function of k for different squared BV frequency N^. The blue 
line represents the N^/£l 2 = — O.f case, and the red line, N^/Vl 2 = — O.Of. The Elsasser number is A = 0.1 and the Roberts 
number is q = 10 -6 . 

place. Note also that the growth rate of the latter is effectively independent of wavenumber for a broad 
range of k. 

In Figs 2 and 3, we present the behaviour of the growth rates as A and iVjj vary. In both case, the 
smaller A and N^, the smaller the double-diffusive growth rate. Numerically, it appears to scale as 




We show this analytically in Section 3.2 with an asymptotic analysis. Being the product of two small 
quantities, realistic growth rates will be much smaller than the orbital frequency, though still dynamically 
significant. For example, taking A = 10~ 3 and N^/fl = 1CP 2 furnishes a growth rate of s ~ 10 _5 fi, and 
thus an c-folding time of some 15,000 orbits. 

The range of unstable wavenumbers associated with the double-diffusive instability is governed by the 
Roberts number q and N^. It extends from, roughly, the resistive scale k n ~ A down to k c ~ lA^/f^Ag -1 / 2 
(in units of fc^). 

3.3. Asymptotic analysis 

In this section we derive an analytic expression for the growth rates in the double-diffusive regime: when 
wavelengths are much less than the resistive scale l v but much longer than the cutoff scale l c . These modes 
are thus heavily influenced by electrical diffusion but are relatively insensitive to radiative diffusion. An 
asymptotic analysis of very long wavelength modes, which connects this paper to the local stability results 
of Balbus (1995), is presented in Appendix A. 

In order to calculate modes in the double-diffusive regime, we set set their lengthscale to a value of order 
the Alfven length, i.e. k ~ fc^. In addition, a small scaling parameter e is introduced so that < e « 1. 
From the outset we assume that both A and N^/Cl 2 are equally small and of order e. Subsequently, the 
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growth rate s is expanded in powers of e, substituted into the dispersion relation, and then terms of equal 
order are collected. This procedure permits the calculation of the leading order terms of all five modes. 

We find that the two hydrodynamic convective modes oscillate at the epicyclic frequency at leading 
(zeroth) order and decay at a rate — A Q. There exist two other modes, on the other hand, which both decay 
at the rate -A -1 (fc/fc^) 2 0. 

The fifth mode is the double-diffusive mode and it possesses a growth rate that scales like e 2 . It can be 
determined by balancing the last two terms in the dispersion relation Eq. (14). We find, to leading order, 

N 2 

s= f AO. (27) 

K 

Stability of the mode is governed by the classical Schwarzchild criterion. So when < 0, the mode grows. 
Note also that the growth rate is independent of the wavenumber k in this regime. The expression (27) 
is, in fact, quite general and holds throughout most of the double-diffusive range of k, as confirmed by 
the numerical solution (Figs 2 and 3). When k becomes small, however, the analysis breaks down and the 
growth rate connects to the classical-MRI curve. And when k is very large, thermal diffusion stabilises the 
mode. 



3-4- Magneto strophic balances 

In this subsection we characterise the instability in a physically illuminating way, by showing how it 
relics on certain steady balances. While the resistive MRI (or magnetostrophic MM) works by balancing 
the Lorentz and incrtial forces in the momentum equation (as discussed in Petitdemange et al. 2008), the 
double-diffusive instability requires in addition the balancing of advection and diffusion in the induction 
equation. Recognising this fact helps simplify some of the algebra of the previous section, and adds support 
to the previous asymptotic estimate. 

3.4-1- Balancing the momentum equation 

Consider the linearised momentum equation (9). Assume the left hand side is subdominant. This means 
that to leading order we have a steady balance between the Lorentz and inertial forces (magnetostrophic 
balance), also relevant to the Earth's core (Petitdemange et al. 2008). By keeping the growth rate explicit in 
the induction and entropy equation, but setting £ = 0, we may derive a (reduced) cubic dispersion relation, 

k 2 s 3 + 2k 2 V k 2 s 2 + \{N 2 R + n 2 )v\ k 2 + (4 + K 2 V 2 )k 4 ] s + Nft nv A k 4 = 0, (28) 

with ft 2 = (dfl 2 / din R) . While an approximation, this equation captures the physical essence of the two 
MRI modes and the double-diffusive mode. We may obtain an analytic expression by solving for k 2 , which, 
after inversion, yields the growth rates of these modes. For purely real s: 



2^ g 2 + (JV 2 + n 2 )v 2^ s VA syjAK^Ws + {N 2 R + Q 2 ) 2 v\ - 4k 2 v\ s 2 
±_ 2{v\ s + nN 2 R v 2 A + K 2 n 2 s) 2(v 4 a s + t 1 N 2 v 2 a + k 2 t 1 2 s) ' ( ' 

There are thus two branches, which we plot in Fig. 4 for representative parameters. The solid blue line 
indicates the 'minus' branch and the dashed red line the 'positive' branch. In addition, we plot the full 
solution to the fifth-order dispersion relation (14) (the dotted black line). As is clear, the two branches of the 
reduced magnetostrophic solution present a reasonable approximation to the actual MRI and double-diffusive 
modes for all k. This strengthens the claim that magnetostrophic balance is central to the mechanism of 
instability in both cases. 

3-4-2. Balancing the induction equation 

It is also instructive to consider the limit in which the time derivatives are negigiblc in the induction 
equation as well. We thus set the left-side of the linearised equation (10) to zero. This means that in 
addition to the steady balance between the Lorentz force and the inertial forces in the momentum equation, 
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Figure 4: Growth rate of the purely real unstable mode s as a function of k as determined from the inverse of (29) . The red 
dashed line indicates the plus branch and the solid blue line the minus branch in (29). The dotted black curve represents the 
full solution to the fifth-order dispersion relation (14). The parameters chosen are N^/Q 2 = —0.1, A = 0.1, fl 2 /Q 2 = —3, 
n/Sl = l. 

there is a steady balance between magnetic diffusion and magnetic advection/distortion. This leaves only 
one time-derivative in the linearised equations (9)— (11). As a result, the growth rate is straightforward to 
calculate: 

-N% T) v\ k 2 



fl 2 v 2 A + (v\ + r) 2 n 2 ) k 2 



(30) 



This expression yields the (scaled) asymptotic growth rate of the double-diffusive instability (27) in the 
limit of large resistivity and for k > va/Q- On the other hand, for sufficiently small k we obtain the growth 
rate of a decaying purely resistive mode described in the Appendix by Eq. (A. 4). These two modes, however, 
occur on different solution branches and the transition between them, in the above, is discontinuous. 

The two steady balances upon which the double-diffusive mode exploits, may be understood as a can- 
celling of the forces associated with the magnetic field and differential rotation. Differential rotation will 
attempt to stretch out a field line, and magnetic diffusion will relax the ensuing magnetic tension. In 
the meantime, however, what magnetic tension is generated can counterbalance the inertial forces in the 
momentum equation. This is important because these inertial forces would otherwise stabilise the mode. 

3.4-3. Angular momentum and heat transport 

The solution obtained in the previous subsection presents a convenient framework to compute the angular 
momentum and heat fluxes aroused by a single linear mode. These may offer an insight into the direction 
of these fluxes in the saturated nonlinear state. 

The thermal heat flux associated with a single Fourier mode of the double-diffusive instability we define 

as 

Jt = 2Re(u' x T') = -2Re(u' x p'), (31) 



where the ovcrline indicates the complex conjugate. From the linearised equations, with s negligible in the 
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momentum and induction equations, we obtain 



-pd R S 
7s 



(32) 



With s > and positive, from Eq. (30), and a negative radial entropy gradient, the double-diffusive mode 
always transports heat outwards, at a rate proportional to the gradient itself. This is as expected: the mode 
will act against the destabilising condition from which it arose. 
The angular momentum flux we define as 



M<K)--r—MKK) 



47T p 



(33) 



Using the eigenfunction of the double-diffusive instability, this can be reexpressed as 



K 



4A 



O 2 



v\ k 2 



(34) 



The double-diffusive mode works on scales k > va/&, and we take A to be small. Thus the first term 
(associated with the Reynolds stress) in the square brackets will dominate the second, which ensures that 
angular momentum is transported inwards. However, for very small k, i.e. for the decaying resistive mode, 
the situation is reversed. 



3.5. Physical description of the unstable mode 

Figure 5 presents a schematic diagram of how the double-diffusive instability works in a physically 
intuitive way, developing the ideas put forward at the beginning of Section 2. Panel (a) shows two fluid 
blobs at different vertical locations but at the same radial location (the violet circles). They are connected 
by a vertical magnetic field line. The two blobs are embedded in a negative radial entropy gradient, but 
they possess the same entropy So as each other, because they inhabit the same radius. Fluid at smaller radii 
(to their left) possesses greater entropy, S+i and S+2, and fluid at greater radii (to their right) possesses 
lesser entropy, S_i and S_2- In panel (b) the blobs are perturbed horizontally, initially deforming their 
magnetic tether. Because of the large resistivity, the blobs will slip through the magnetic field, and so the 
magnetorotational instability does not set in. Nevertheless, the field will permit some amount of angular 
momentum to be exchanged between the two blobs and they adopt new orbits at their new radii. But the 
displaced fluid blobs now possess a different entropy to their surroundings: too little entropy in the case of 
the upper blob (So < S+i), and too much in the case of the lower (So > S_i). Unless thermal diffusion acts 
rapidly to equilibriate them with their ambient gas, they cannot sustain their current positions. Instead they 
will continue moving radially (pink arrows). This motion will bring them into contact with new diffusing 
magnetic field lines, panel (c), which will tether them to new fluid partners (the grey circles). Now the new 
magnetic field lines will facilitate another transient torque, which will transfer more angular momentum. 
As earlier, the violet blobs adopt new radii but the disparity between their entropy and their environments' 
entropy is even larger than before, panel (d), and they are compelled to continue drifting. The cycle 
continues, leading to instability. Notice that the cycle incorporates the B field, resistivity, and dS/dR < 0, 
three ingredients which are represented in the stability discriminant in ao, Eq. (19). 

The essential features of the instability are: the suppression of stabilising rotation by the diffusing 
magnetic field, and slow (or negligible) thermal diffusion. This justifies calling the instability double-diffusive: 
angular momentum is diffused rapidly, eliminating the stable angular momentum stratification. But heat is 
diffused slowly, allowing the unstable entropy stratification to do its work unhindered. 

ft can be shown that a sufficiently large viscosity can do the same job as the magnetic field here: and, 
in fact, a hydrodynamical disk in which -C 1 exhibits an analogous instability (where v is kinematic 
viscosity). Consequently, we may interpret the instability as the 'dual' of the Goldreich and Schubert 
instability (Goldreich and Schubert 1967), which relies on the negation of a stabilising entropy gradient, 
by rapid thermal diffusion, and the tapping of an unstable angular momentum gradient (unimpeded by the 
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Figure 5: Four panels showing schematically the physical mechanism of the double-diffusive instability. Panel (a) presents 
two fluid blobs, embedded in a negative entropy gradient: S+2 > S+i > So > S—i > S— 2- The fluid elements are separated 
vertically but connected by a magnetic tether. The following panels, in order, show the development of a horizontal perturbation 
in time, under the influence of speedy angular momentum transport (by the diffusing magnetic field) and the radial entropy 
gradient. See Section 3.5 
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weak viscosity) . The dead-zone double-diffusion instability is also related to salt- fingering in the ocean: and, 
in fact, the linear stage of a single unstable mode will be characterised by a vertical sequence of horizontal 
'fingers' or channel flows, much as in the MRI. However, these flows, unlike classical MRI, are not nonlinear 
solutions, and will quickly break down into some form of disordered mixing. Another important difference 
is that nearly all the unstable modes possess comparable growth rates, so no one mode will dominate. 
Noticeable channel flows, or fingers, are therefore unlikely to be observed in simulations. 

4. Conclusion 

This paper introduces a novel linear instability which occurs on short length-scales and which could be 
at work within the dead-zones of protostcllar disks. The instability is double-diffusive and relies on (a) 
the energy provided by a weak negative entropy gradient, (b) a diffusing magnetic field to relax angular 
momentum conservation, and (c) very slow thermal diffusion. It grows on a range of lengthscales extending 
from the resistive length down to a length related to the thermal scale. Across this potentially large range 
the growth rate scales like s <~ A (N^/fl 2 ) ft. Though the unstable mode will almost always physically fit 
into the dead-zone (unlike the MRI), it may grow quite slowly. Nevertheless it may give rise to dynamically 
interesting velocity amplitudes after some 10 4 orbits (depending on the parameters). In principle, the 
instability's nonlinear saturation may play a significant role in the dynamics of dead-zones, leading, in 
particular, to radial and vertical mixing, which in turn may influence dust grain dynamics (settling and 
clumping) on one hand, and the global temperature structure of the disk, on the other. 

However, a number of issues need to be settled, and these form the basis for future work. First, the 
governing parameters of the system need to be better constrained, so that one can be assured the instability 
grows on dynamically meaningful timescales. This requires realistic estimates for the magnetic field strength, 
the radial entropy gradient, and the resistivity. Second, the linear analysis should be extended to the peculiar 
nonideal MHD conditions present in the interior of protostellar disks, where the instability will be altered by 
Hall effects and ionisation and dust grain chemistry. Third, the interaction between the instability and much 
faster processes in the active zones needs to be clarified. For instance, the turbulent motions in these regions 
may transport the large-scale magnetic field on a similar time to the growth time (Fromang and Stone 2009, 
Lesur and Longaretti 2009). On the other hand, rapid density waves, stimulated by the turbulence, may 
also interact with the instability (Fleming and Stone 2003, Heinemann and Papaloizou 2009a, b). Fourth, 
the instability's nonlinear saturation should be simulated numerically in order to establish whether it leads 
to sustained disordered flows. 

Finally, this instability is not restricted to protostellar disks, as the necessary ingredients for instability 
are rather generic. In particular, double-diffusive modes may be excited within the fluid cores of rapidly 
rotating planets and protoplanets, possibly modulating pre-existing dynamo activity, or perhaps directly 
instigating convection itself in bodies that spin so rapidly that regular buoyancy instabilities are quenched. 
These issues we address in future work. 
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Appendix A. Asymptotic description of long wavelength modes 

In this appendix we calculate the growth rates of the five modes of Eq. (14) in the limit of long wavelength, 
i.e. on scales longer than the resistive scale l v . As earlier, we adopt O -1 for the unit of time, and the Alfvcn 
wavenumber (Q/va) for the unit of wavenumber. In addition, we drop the thermal diffusion terms, assuming 
that k and £ are so small that these effects will be negligible. The dispersion relation (14) in these units is 

s 5 + b A s 4 + 6 3 s 3 + b 2 s 2 + 61 s + b = 0, (A.l) 

with coefficients 

b 4 = 2A- 1 fc 2 , 

b 3 = k 2 + N^ + 2k 2 + A' 2 k 4 , 

b 2 = 2A- 1 k 2 (K 2 + N 2 +k 2 ), 

h = k 2 (k 2 - 4) + fc 2 (l + A" 2 fc 2 )(A^ 2 + k 2 ), 

fe = A- 1 iV 2 fc 4 , 

and in which both k 2 and have been scaled by O 2 . 

A small scaling parameter e is introduced so that < e <C 1. From the outset we assume that both A 
and are equally small, and let A ~ ~ e. In order to pick out the longest wavelength modes we set 
k <~ e 2 , which means that we examine lengthscales A -1 times the resistive scale. Next the growth rate s is 
expanded in small e, which allows the computation of each mode analytically to leading order 

The two hydrodynamical convective modes possess an expansion like s = s + Si e + ,s 2 e 2 + • ■ • , with the 
leading order term given by s — ±i k. Thus the modes oscillate at the epicyclic frequency at leading order. 
The real part of the growth rate comes in only at seventh order in e 

Ak 4 

Re( S )---^- + 0( £ 8 ), (A.2) 

and thus decays at a rate proportional to 77. Note that if a different scaling had been adopted in which Aj. ~ 1 
then the leading order term would be so = ±y / —n 2 — iV| and stability would follow the Hoiland-Solberg 
criterion, k 2 + > 0. 

The leading order term of the two MRI (or convective-MRI) modes scales as e 2 and is determined by 
balancing the second and fourth term in Eq. (A.l). We find 

»=4fI5+ o <« 4 >< < a - 3 » 
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which returns the MRI stability condition d£l 2 /dlni? > 0. Thus these very long scales will be MM unstable. 
Note that the entropy gradient has not entered at leading order, it being too small to determine the stability 
properties of the gas. But if we had assumed that <~ 1 the growth rate to leading order would be 
proportional to y 7 — k 2 — iVj| + 4, and the instability criterion would reproduce that for convective-MRI 
modes (Balbus 1995, see also Papaloizou and Szuszkiewicz 1992). 

The last mode is a simple resistive mode whose growth rate scales like e and which can be determined 
from the last two terms of the dispersion relation. We find 

k 2 N 2 

— Twmm +0{e ° ) - lAA> 

and thus the mode decays at a rate proportional to 77. Note that this mode is not the double-diffusive mode 
computed in the main body of the paper. It exists on a separate branch of the dispersion relation and will 
always decay for all k. However, just like the double-diffusive instability, it establishes magnetostrophic 
balances in both the induction and momentum equation (see Section 3.4.2). 
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